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Abstract 

With the advent of high precision neutrino scattering experiments comes the need for improved 
radiative corrections. We present a phenomenological analysis of some contributions to the produc- 
tion of photons in neutrino neutral current scattering that are relevant to experiments subsuming 
the 1% level. 
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I. INTRODUCTION 



Recent neutrino scattering experiments, particularly those searching for electron neutrino 
appearance [H El E] by oscillation from muon neutrino beams, report signals below the 1% 
level. Such unprecedented levels of sensitivity demand corresponding efforts to determine 
backgrounds. Radiative corrections are clearly expected at this level. 

Many years ago, phenomenological analyses of radiative corrections were applied to stud- 
ies of pion decay [H [5], to constrain corrections to the weak interaction from non-V-A 
currents. Their accuracy was confirmed by precise calculations when the Standard Model 
became established. 

Although many contributions are possible, obviating unambiguous predictions, similar 
phenomenological analyses can provide overall scale, energy and angular distributions for 
photon production in neutrino scattering. Such methods were used in Ej to help explain 
the Gargamelle single photon signal [S] . Armed with these predictions, experimental results 
can then be used to determine whether the total of such contributions are significant within 
each experiment. 

Of course, Monte Carlo calculations by experimental groups include bremsstrahlung from 
charged particles in the final state and radiative decays of N* and A resonances. However, 
these are s-channel effects. We consider here photons produced from the particles being 
exchanged in the t-channel, and so less obviously connected to the external line quanta. In 
general, such contributions are necessary simply to maintain gauge invariance. 

We provide here two examples of amplitudes relevant to neutral current neutrino scatter- 
ing, which are especially important when distinguishing final state electrons and photons is 
difficult. We eschew details such as interferences between different contributions, as the final 
state distributions in our examples differ significantly, minimizing the overlap. (In addition, 
it is usually the case that possible interfering contributions involve the exchange of much 
higher mass particles with much smaller phenomenological couplings [9]; hence, interference 
effects are expected to be minimal.) Our goal is simply to provide experimental groups with 
distribution formulas that they can use to extract these contributions or determine limits 
on them. 



We provide results for both neutrino and antineutrino scattering. Our results may be 
easily extended to other similar contributions both in neutral and charged current scat- 
tering. Although we focus on modest energy studies, our results are fully relativistically 
covariant and thus may be applied to any energy. Additional contributions or Regge trajec- 
tory generalizations of the meson exchanges are likely necessary for application to very high 
energies. Our results are also straightforwardly generalizable to other related processes such 
as charged current interactions, which we will present elsewhere. 



II. LAGRANGIAN 



Photon production in neutral current neutrino scattering produces a background to the 
identification of electron neutrino events in (generally non-magnetic) detectors which distin- 
guish poorly between electrons and photons. Standard radiative corrections, such as those 
due to bremsstrahlung of a photon from charged particles in the target or final state, have 
been closely examined previously [10]. See also [Tij IT2] and references therein. Additional 
processes, such as production of a A baryon or N* followed by its decay back to a nucleon 
and a photon [13], are also already included in many experimental Monte Carlo assessments 
of backgrounds [9] . 

Recently, however, a "new" triangle anomaly has been identified in reference |14j . here- 
inafter referred to as H3, which can contribute a previously unconsidered source of such 
background. However, the structure of that contribution is very similar to a phenomenolog- 
ical contribution which we display in Fig. ([!]) . Here the Z-boson carrying the neutral current 
interaction from the neutrino line mixes into a vector boson, in the familiar fashion of Vector 
Meson Dominance (VMD). [15] In this case, the hadronic vector meson is an uj (or a p°) 
which undergoes a virtual decay to the photon of interest and a pion in the t-channel. This 
last couples strongly to the hadron target (nucleon or nucleus) completing the interaction. 

The advantage over the H3 approach is that the vertex strength is known phenomeno- 
logically from the rate of the decay 

UJ — > 7T° + 7 

and similarly for the p° case. The pion coupling to the nucleon is also well-known [T5]. The 
strength of the Z — uj mixing is determined by VMD and the off-shell variation of this mixing 



is easily determined by a well-known extension [TTl [18] , originally worked out for the case of 
isospin violating contributions to the nuclear force. We explicitly calculate the appropriate 



analog in Appendix (A 2) for completeness. 

No other parameters are required, so the prediction of the contribution to the total cross 
section for producing a final state photon is absolute for this graph. We emphasize that the 
structure of the central vertex has the same vector-vector-axial- vector coupling structure as 
in H3, due to the axial vector nature of the pion current. However, in our approach, all of 
the strong interaction corrections (higher order in QCD and quark-antiquark "resonance" 
effects) relevant to the triangle graph of H3 are fully taken into account phenomenologically. 
Of course, other similar contributions occur with vector- meson recurrences, etc., but these 
predominantly affect only the overall strength, and furthermore, experience has shown that 
at the modest energies of the experiments of interest (LSND [TJ and MiniBooNE [21 E]), the 
sum over all such contributions is likely to be dominated by these leading ones. 

The interaction Lagrangian terms needed in our approach are: 

£/ = eg^e^u^d^F^ + g nNN ^d^ ■ r * (1) 

where \1/ is the nucleon field and f is the vector of usual isospin generators. 

We also make use of the usual Standard Model weak interaction couplings of the Z- 
boson to neutrinos and to quarks. The last are needed to compute the off-shell variation 
of the VMD mixing between the Z and the u. The uj — 7r — 7 coupling constant, g W7r7 = 
eg u (with the explicit factor of the electromagnetic coupling, e, stripped out for clarity) is 
determined from the experimental value of the (on-shell) radiative uj — > it + 7 decay rate, 



shown in Appendix (Al). We will discuss possible off-shell effects on this compared with 



the corresponding effects for the H3 approach in our discussion section (IV). 



III. CROSS SECTION 



We consider the process shown in FigjTJ In terms of the labeled momentum 4-vectors, 
the squared amplitude is 

~ (g 2 -M!) 2 (g 2 -M2)2(Q2_ M 2 )2 ^ W Pf M N ) [{k t t ) + {k t f ) ) , {2) 
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FIG. 1: Specific diagram considered in this analysis. Variants are discussed in the text. 

where the neutrino- Z coupling is given by g v z = (<?/2) cosOw in terms of the SU{2) gauge 
coupling g. The remaining coupling constants will be discussed in what follows. In the 
center of mass (CM) frame, the momenta can be written explicitly as 



ti = (Ee i: E (i ) (3) 

Pi = (E Pi ,-E ei ) (4) 

I, = (E t/ ,E tf ) (5) 

Pf = ( E p r P P f) ( 6 ) 

k = (E k ,E k ), (7) 



where we employ the shorthand E { = pi to indicate a massless particle's 3-momentum. 
To obtain the desired cross-section, we first partially evaluate the phase space integrals by 
making use of momentum conserving delta functions. We find the differential cross-section 
for the photon in the final state is 

dcr M%gl z gl^gl NN r gl z {q 2 )q 2 Q 2 ({k ■ ttf + (k ■ i f ) 2 ) 

dE k dfi {2^YE ti {E ti + E pi ) J ^{q 2 - M z )\q 2 - M 2 ) 2 (Q 2 - 1 ' 

in terms of its energy (E k ) and opening angle from the beam direction (// = cos#). Here the 
t-channel momentum transfers are: q 2 = — 2E ii Ei f (l — /i£ f ) and Q 2 = q 2 — 2k ■ t; L + Ik • if, 



where fii f is the cosine of the opening angle between the neutrino in the final and initial 
state. It is related to /x and the cosine \ii jk of the opening angle between the photon and 
the final state neutrino by 



fl tf =/!/!^fc + yl-/! 2 yl-/if /fc COS0. (9) 

This is the only source of <fi dependence in the system. Momentum conservation fixes 

fii f k = (v^ 2 - 2yfi(Ee f + E k ) + 2E k E if - M 2 N ) (10) 
2E if E k v 

as a function of the invariant y/s = E £i + E Pi . Requiring that yn fk G { — 1, 1} yields the Ee f 
integration limits 

in y/s 2 — Mji — 2^SE K 

E '' = 2-fs (U) 

y/s 2 — Mf] — 2^sE K 

Ee ' = 2V^lE k • (12) 

Additionally, we require E k < ^/%{yfs 2 — M^) to maintain positive values for these energies. 
The uj — Z mixing coupling function, g 2 z (q 2 ), is found from the self energy diagram in Fig|7] 



of Appendix (A2). It is almost constant in the space-like regime where it is needed here. 
Taking the limit \q 2 \ <C M§ and \Q 2 \ ^> M 2 , the cross section simplifies to 

^_ 2 2 2 2 2 

da MNQuZdiuZOuj-yndnNN { AT? //, « \2 , // o \2\ f i- Q l-x n\ 

md, = { 2«y El ,(E tl + ejmi J dE " «* ■ ll) + <* ' «'> ) 7 (13) 

where we have neglected the g 2 dependence of and assume an average value. Integrating 
over 6, we obtain 



da M 2 N E 2 gl z gl z gl^gl NN 



J dE lf (E\{\ - fi) 2 + E 2 f (l - fi ifk ) 2 ) (14) 



dE k d(i {2^fE ti {E h +E Pi )M z 

1 { a — b c 3 — 2ac 2 + abc — b + a 



/ 2 (&-c) 2 \(6 2 -l)l (c 2 -l) 
where we have introduced 



/ = 2E ii E^\-p? s J\- V ? tfk (15) 
a = J — — (16) 



2E ti Ee f {l - nni f k) + 2E k E u {l - /i) - 2E k E £f {l - fi £fk ) 

f 

2Ee t E ef (l - ^ £fk ) + M 2 
f 



(17) 
(18) 



200 MeV 




FIG. 2: Lab frame cross section with beam energy Eg. = 200 MeV. 



to simplify the notation. The dimensionless quantities a, b and c are all greater than unity, 
assuming physical parameter values. This leaves only the one- dimensional Eg, integral to 
perform. 



A. Cross Section Characteristics 

In the CM form, the predicted cross section is weakly peaked in the backward direction 
with an energy maximum near the highest allowed energies due to the overall factor of E^. 
Boosting to the lab frame pushes the angular distribution forward and spreads the energy 
of the photon in the standard way. Numerically integrating Eq. Q, we plot the differential 
cross-section for beam energies of 200 MeV, 350 MeV, 500 MeV and 1000 MeV in Figsj2| 
|3j [4] and |5j respectively. These are boosted to the lab frame for convenience. In each case, 
we display the E n and cos 9 dependent contour plots as well as energy projection panels 
for the result of integrating over one of the variables. The total cross section is also noted 
for reference. The angular distribution moves toward the forward peak with increasing 
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FIG. 3: Lab frame cross section with beam energy Eg. = 350 MeV. 
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FIG. 4: Lab frame cross section with beam energy E^ = 500 MeV. 




FIG. 5: Lab frame cross section with beam energy Eg. = 1000 MeV. 

neutrino energy due to the growing boosts from the CM to the lab frame. The distribution 
consistently peaks near the center of the kinematically allowed photon energy range. 

Integrating over the final state photon energy and angular distribution, we plot the total 
cross section as a function of neutrino beam energy in Fig|6] At high energies (E u ^> Mn) the 
cross section grows as \[E~ V . Near threshold, it grows as E%, as can be seen in the logarithmic 
low energy insert plot. In the energy region below 1000 MeV, the cross section does not 
exceed 10~ 3 fbarns, which is roughly three orders of magnitude less than the typical charged 
current cross sections at these energies [19] . Nevertheless, it is still at a level where it can 
affect current 0, |3] and next generation |20j experiments that probe neutrino interactions 
with sub-percent sensitivities. Current and proposed long baseline oscillation experiments 
(see, for example pTl [221 [231 [21] and references therein) will be sensitive to this class of 
processes with the order of magnitude enhanced cross sections shown in the high energy 
region of FigJHj 




FIG. 6: Total cross section as a function of neutrino beam energy 
IV. DISCUSSION 

Other processes, related to that in FigjTJ such as by exchange of the uj and 7r° for other 
mesons (with the same quantum numbers), will contribute to similar production of photons 
in the final state. These will differ from our calculation here only by the meson masses and 
coupling constants. 

We also observe that this process yields identical results for both neutrino and anti- 
neutrino scattering. The only difference between these resides in a sign change between the 
axial-vector neutrino-Z coupling. The asymmetry term vanishes when Lorentz contracted 
with the rest of the diagram, which is symmetric under the free indices of the vertex. (We 
have checked this feature by direct calculation.) This is not the case for non-radiative neutral 
current neutrino scattering, where there is a difference between neutrino and anti-neutrino 
amplitudes. Additionally, we find by means of direct computation that many variants of 
Figjl] vanish due to similar symmetry reasons. In particular, 



Amplitudes from diagrams with axial vector mesons replacing the vector u vanish. 



• The amplitude of the process with the ir and u mesons interchanged is zero. 

• Amplitudes from "reversed" diagrams (where the 7r couples to the neutrino) vanish. 

We have used a phenomenological on-shell coupling strength for the u — tt — 7 vertex, 
which is a commonly used phenomenological approach. The slow off-shell variation of the 
computed Z — u mixing supports such an approach, although a three-particle vertex may 
well behave differently. The variation one expects is for a decrease in amplitude as vertex 
form factors act to suppress the effective coupling [25J . This may be ultimately recovered by 
summing over higher mass Regge recurrences in an extension of our approach, which would 
eventually transition at very high energies to a purely quark picture of the couplings. 

Such a view raises the question as to whether the H3 axial anomaly might not give a 
more reliable evaluation of the strength of these processes after all. However, the use of 
the anomaly is only guaranteed accurate when currents coupling to the vertices are point- 
like. In the case of the pion, the axial current satisfies this condition, but for the u, being 
a composite object, the vector isoscalar current cannot provide the same guarantee. The 
vertex structure for uj — q — q coupling produces the same kind of uncertainties in the H3 
approach as we encounter here; perhaps more so as the u momentum wavefunction in the 
quark basis may only be known via Lattice QCD calculations, although crude quark model 
calculations do produce a 3-ir hadronic width consistent with experiment. [26J We believe 
the issue is less serious for our calculation as the Z — uo mixing is effectively used only within 
a region of order only a few mass squared units from the on-shell point. 

Our view is that for low to modest energy scattering, such as is relevant to the experiments 
discussed here, our phenomenological approach is as valid as any other, although it is most 
reliable for momentum transfer and angular distributions, and least reliable for absolute 
rates. The most certain point is that these t-channel processes will be the same in neutrino 
and in anti-neutrino scattering, so a comparison of the two running modes for MiniBooNE 
will set the most stringent limits on the strength of these contributions to the total rate. 

Another significant background also occurs when a neutral pion is produced instead of a 
photon, but one of the photons from the decay of this pion is lost to the detector. Fortunately, 
in this case the event rate can be normalized to the corresponding case in charged current 
neutrino scattering, which produces a neutral pion in conjunction with the charged lepton. 



Although this may dominantly occur due to intermediate state processes, such as production 
of a A baryon or N* followed by its decay back to a nucleon and a pion (see [27] and reference 
therein), concern also arises regarding other processes, including those that may be coherent 
over the entire nuclear target with attendant amplification of the rate [281 129] . The approach 
taken above can also be applied to such calculations where again, the pion is not produced in 
an s-channel fashion from a nucleon or one of its N* or A excited states, but from a t-channel 
process, corresponding to the photon production above. We will discuss that process in a 
separate paper. 

Finally, we comment on the effect of replacing the u by a p°: For the case of p meson 
exchange (in place of the a>), the analogous plots to those displayed for u;-exchange are 
very similar in shape and absolute magnitude. The ratio of separate contributions depends 
upon both the ratios of the vector meson-pion-photon coupling constants, g'y_ 7 _ 7r , estimated 
in the Appendix, as well as the vector meson-Z coupling constants, gvz- We estimate 
a p /er w m ((7 p _ 7 _ 7r /(7 a ;- 7 - 7r ) 2 x (g p z/gujz) 2 ~ 1-6, using parameter values taken from Ref . ( [30] ) . 
The suppression arises primarily from the coupling associated with the meson— it — 7 vertex 
while the enhancement is due to Z-meson mixing. See the Appendix [A] for more details. 

The amplitudes for these processes are comparable, and significant interference is ex- 
pected to be possible. From this effect the overall cross sections may be modified by a factor 
between 0.07 and 5.1 for total destructive and constructive interference respectively. The 
cross section yielded by the lower limit is well below current or expected future experimental 
sensitivities and therefore forms a negligible background in that case. These considerations 
are interesting for future work, because of the parallel process of coherent pion production: 
If the interference is destructive, there can be a significant difference between coherent pion 
production between charged and neutral current this interference cannot occur in 

the charged current case. 

In conclusion, we remark that whether the interference is constructive or destructive, our 
results apply equally to neutrino and to antineutrino neutral current scattering, since that 
only affects the overall sign of the amplitude by the neutrino coupling to the Z-boson. Thus, 
the process discussed here must contribute equally to neutrino and to antineutrino neutral 
current scattering, and with the same energy and angular dependence for the appearance of 
the photon. 
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APPENDIX A: SUBSIDIARY ISSUES 

For completeness, we show the explicit calculations of the decay rate for the uj and the 
Z — uj mixing amplitude. 



1. T{uJ — > 7T + 7) 

We calculate the uj decay width r(u; — > n + 7) from the phenomenological vertex Cj = 
eg^^e^^uj^d^Tr F^ u . This same form arises in the triangle anomaly considered in H3. 
Neglecting the ir° mass, the squared amplitude of this process is 

A 2 = _^TL k . q= e Wg* (Al) 
3 6 

where k and q are the photon and pion momenta, respectively. Evaluating the decay width 
yields 

aq 2 M 3 

r(u,^7r + 7) = -?^. (A2) 



Fitting this to the observed decay width [50], we extract the coupling constant 
9uyn — 1.8/M^, which we use in the numerical examples throughout our analysis. Evaluat- 




h-k 



FIG. 7: Self energy diagram leading to Z — uo mixing. 



ing the tree diagram without the Lorentz structure of the anomaly yields g W77r ~ 1.2 /M^ 



i.e., a factor of 2/3 smaller than is obtained from Eq. (A2). These are written in the form of 
a dimensionless coupling relative to the uo mass to help illustrate the scale associated with 
the effective vertex and strong coupling strength. A similar exercise may be performed with 
the p° decay, in which case one extracts <? P77r = 0.55/M p . 



2. p 2 dependence of the Z — uo mixing amplitude 

The p 2 dependence of the Z — uo mixing coefficient g w z{p 2 ) may be found from the self- 
energy diagram shown in Fig |7j Following references ( [I7J [18] ) , we parameterize the oo — qq 
vertex by the form factor g uqq -M 2 /(M 2 — p 2 ), where the g uq g ~ 3.1 coupling is extracted 
from uo — > 7T° 7r + 7r~ decay measurements [30], under the usual quark model assumption 
that inclusive processes may be well-approximated with the use of quark-hadron duality. 
The scale, M, defines the finite substructure of the uo meson. Evaluating this diagram, we 
find the self-energy 

IV = ^rp^r dW I dxdydz5(x + y + z - 1) (A3) 
3(27r) d cos6'vi/ Jo 

d I- 3 , 

{t 2 — p 2 z(z — 1) — m 2 (y + z) — M 2 x^j 

where the d-dimensional integral is written to aid in the renormalization of the divergent £ 2 
term via dimensional regularization and m q is the average light quark mass. Considering only 
the terms that contribute to the p 2 dependence of uo — Z mixing and dropping logarithmic 
contributions, we find 

2 -gg^M 2 sin 2 6 W yi r^-* p 2 z(z-l)+m 2 q 
g.zKV)- 127l 2 cos9w J Q az J %2^_ 1)+m 2 + x(M 2_ m 2)- ^ 



Taking the reasonable limit m 2 <C M 2 , p 2 , we evaluate the above integral to yield 



2 ^ gg ,M 2 sin 2 fl w 1_ f _ + _ ^ ^ ^ / j_jA 

727H cos tV 77^ [ \ ~V J 

+ e ^7(4 - 77) + (4 - 377) ln(l - 77) - 77 2 (6 + 77) In (^^j ) + ( e2 )| > 

with the dimensionless quantities, 77 = p 2 /M 2 and e = m 2 /M 2 . For the p 2 < range 
relevant to this analysis, g w z{p 2 ) is real. The e correction term is well-behaved and negligible 



for realistic parameter values where e ~ 10~ 5 . The leading order term of Eq. (A5) breaks 
down for — p 2 < 1 MeV 2 or —77 < 1.6 x 10~ 6 . The limiting behavior of g^zip 2 ) is easy to 
extract and is useful in estimating the physical parameter range 

+ +0 (J) (-„»!) 

24n 2 cos6 w I A(l — 77ln(— 77)) + 0(r]) (-77 < 1) 

For the neutrino beam energies studied here between 200 — 1000 MeV, we find that the 
average coupling g^z = 600 MeV 2 , under the assumptions that m q ~ 3 MeV [30J and 

M ~ M,„. 



3. Relative strengths of Z — p and Z — lo mixing 

On the basis of SU(3) flavor symmetry, one expects a similar result for the p°, except for 
the effect of isospin. In Fig. ([7]), the u> couples equally to the u and d quarks that contribute 
to the loop, which weights the Zqq couplings, viz. 

9zuu = . 9 a (^sin 2 ^-l) (A7) 
4 cos &w o 

9zm = ^^(l-^sin 2 M (A8) 
4 cos &w o 

equally, whereas the p° does so with opposite signs due to isospin. (The 75 parts of the 
couplings are omitted as irrelevant; they cannot contribute to the mixing since there is only 
one momentum available to combine with the 4-dimensional Levi-Civita tensor that comes 
from the fermion trace.) As might be expected from the fact that the Z-boson is dominantly 
isospin 1 like the p, the Z — p mixing is enhanced relative to the Z — u mixing. 



Combining these considerations with the slight phenomenological deviation from SU (3) 
flavor symmetry yields 

9pZ _ 9zdd ~ 9Zuu x Qpqq (A9) 
9uZ gzdd + 9Zuu Qwqq 



T(p — > 7T7r) — > 7r7T7r) 



1-2 sin 2 0^ 
2 sin 2 9 W ^ r(c<j — > nnn) 0(p — > 7T7t) 
= 4.1, 

where <ft denotes the phase space integral for the decay, which suppresses the ratio by ap- 
proximately yAn. We take the numerical value of sin 2 6\y ~ 0.2396 from the Particle Data 
Group extrapolation to small four- momentum transfer [30]. The effect of this enhancement 
is included in our final discussion. 



